home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CTRCO.z / CTRCO
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCTTTTRRRRCCCCOOOO((((3333FFFF))))                                                            CCCCTTTTRRRRCCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CTRCO   - CTRCO estimates the condition of a complex triangular matrix.
  10.  
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  13.       SUBROUTINE CTRCO(T,LDT,N,RCOND,Z,JOB)
  14.  
  15. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  16.      On Entry
  17.  
  18.      TTTT COMPLEX(LDT,N)
  19.         T contains the triangular matrix.  The zero
  20.         elements of the matrix are not referenced, and
  21.         the corresponding elements of the array can be
  22.         used to store other information.
  23.  
  24.      LLLLDDDDTTTT INTEGER
  25.         LDT is the leading dimension of the array T.
  26.  
  27.      NNNN INTEGER
  28.         N is the order of the system.
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         = 0         T  is lower triangular.
  32.         = nonzero   T  is upper triangular.  On Return
  33.  
  34.      RRRRCCCCOOOONNNNDDDD REAL
  35.         an estimate of the reciprocal condition of  T .
  36.         For the system  T*X = B , relative perturbations
  37.         in  T  and  B  of size  EPSILON  may cause
  38.         relative perturbations in  X  of size  EPSILON/RCOND .
  39.         If  RCOND  is so small that the logical expression
  40.         1.0 + RCOND .EQ. 1.0
  41.         is true, then  T  may be singular to working
  42.         precision.  In particular,  RCOND  is zero  if
  43.         exact singularity is detected or the estimate
  44.         underflows.
  45.  
  46.      ZZZZ COMPLEX(N)
  47.         a work vector whose contents are usually unimportant.
  48.         If  T  is close to a singular matrix, then  Z  is
  49.         an approximate null vector in the sense that
  50.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  51.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  52.      Subroutines and Functions BLAS CAXPY,CSSCAL,SCASUM Fortran
  53.      ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.